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COSMIC  RAY  EFFECTS  ON  MICROELECTRONICS 
PART  II:  THE  GEOMAGNETIC  CUTOFF  EFFECTS 

1.0  Introduction 


The  ability  of  charged  particle  radiation  to  penetrate  Into  the 
magnetosphere  from  outside  Is  limited  by  the  earth's  magnetic  field. 
Particles  having  low  magnetic  rigidity  (l.e.  momentum  per  unit  charge)  are 
preferentially  turned  back  by  the  field,  so  they  are  unable  to  penetrate 
beyond  some  depth  in  the  magnetosphere.  For  each  point  In  the 
magnetosphere  and  for  each  direction  of  approach  to  that  point,  there 
exists  a  threshold  value  of  magnetic  rigidity,  called  the  geomagnetic 
cutoff.  Below  this  value,  no  charged  particle  can  reach  the  specified 
point  from  the  specified  direction.  Above  this  cutoff  value,  particles 
arrive  at  the  specified  point  from  the  specified  direction  as  though  the 
magnetic  field  were  not  present  at  all  (Lemaitre  and  Vallarta,  1933). 


2.0  The  Geomagnetic  Cutoff 

The  geomagnetic  cutoff  was  first  computed  for  a  perfect  dipole  field  by 
Stormer  (1930).  This  result  Is  given  below: 

P  ■  [l-(  1-cosy  coAW/Cco,,  cosx]2  (1) 

c  / 


for  positively  charged  particles,  where 

Pc  *  magnetic  rigidity  cutoff  in  GeV/ec  (or  GV  for  short), 
r  ■  radial  distance  from  the  offset  dipole  center  In  earth  radii, 
x  ■  geomagnetic  latitude  in  offset  dipole  coordinates, 

and, 

y  ■  the  direction  from  which  the  particle  arrives  with  respect 
to  local  east  In  offset  dipole  coordinates. 


As  we  can  see  from  Stormer's  equation,  Pc  depends  strongly  on  the 
geomagnetic  latitude.  This  means  that  the  geomagnetic  cutoff  will  vary 
drastically  around  the  orbit  of  a  spacecraft,  especially  at  high 
inclinations. 

Stormer's  equation  assumes  the  earth's  field  is  a  perfect  dipole  and 
Ignores  the  presence  of  the  solid  earth.  For  these  reasons,  Stormer's 
equation  does  not  provide  accurate  values  of  the  geomagnetic  cutoff. 

More  accurate  values  of  the  geomagnetic  cutoff  have  been  obtained  by 
Shea  and  Smart  (1975).  These  authors  have  used  an  accurate  model  of  the 
earth's  magnetic  field  and  ray  tracing  techniques,  to  obtain  values  of  the 
vertical  geomagnetic  cutoff.  Their  calculations  take  into  account  the 
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presence  of  the  solid  earth.  Shea  and  Smart  report  the  values  of  the 
vertical  geomagnetic  cutoff  at  20  km  altitude  on  a  world-wide  grid  of 
points.  These  results  are  for  geomagnetical ly  quiet  conditions. 

When  a  large  solar  flare  occurs,  It  usually  produces  a  magnetic  storm 
at  earth.  These  storms  have  a  tendency  to  depress  the  geomagnetic  cutoff 
as  discussed  In  Adams  et  al.  (1081),  who  showed  that  a  typical  suppression 
In  cutoff  could  be  described  by: 

P storm  *  Pet1-0  “  0'54  exp( «PC/2.9GV)]  (2) 

When  the  geomagnetic  cutoff  transmission  function  Is  being  computed  for 
the  purpose  of  modulating  solar  flare  spectra,  eq.  (2)  will  be  used  to 
adjust  the  geomagnetic  cutoff  values  to  storm-time  conditions.  Obviously 
the  actual  degree  of  supresslon  will  vary  from  one  magnetic  storm  to 
another.  The  data  base  at  present  Is  Inadequate  to  characterize  this 
variation. 

As  can  be  seen  from  eq.  (1),  the  geomagnetic  cutoff  varies  with  the 
particle's  arrival  direction  The  vertical  cutoff  (y  -  90  degrees)  Is 
close  to  the  average  value  of  Pc,  averaging  over  all  directions.  Figure 
1  shows  the  results  of  a  detailed  calculation  of  the  geomagnetic  cutoff 
transmission  averaged  over  all  directions  about  a  point  at  20°  latitude, 

0°  longitude  and  400  km  .  The  arrow  Indicates  the  value  of  the  vertical 
cutoff  at  this  location.  The  vertical  cutoff  will  be  used  here  to  apply 
to  all  cosmic  rays  approaching  a  spacecraft  regardless  of  their 
direction.  The  results  of  Shea  and  Smart  (1975)  are  for  20  km  altitude. 

To  estimate  the  cutoff  at  higher  altitudes  we  will  extrapolate  using  eq. 

(1). 


Finally,  as  Adams  et  al .  (1981)  have  discussed,  the  earth  casts  a 
cosmic  ray  shadow  on  a  spacecraft  In  earth  orbit.  The  solid  angle 
of  the  shadow  cast  of  the  solid  earth  Is  approximately, 

Omega  *  2n  {  1-  [(R0  +  h)2-  R2  ]  1/2/  (R0  +  h)f  ,  (3) 

where  Re  Is  the  radius  of  the  earth  and  h  Is  the  spacecraft  altitude. 

This  shadow  Is  cast  on  one  side  of  the  spacecraft  at  any  point  on  the 
orbit.  In  the  results  reported  here,  we  will  neglect  the  earth's 
shadow.  This  will  lead  to  direction  -  averaged  flux  estimates  for  the 
radiation  reaching  the  skin  of  the  spacecraft  that  are  somewhat  high 
(never  more  than  a  factor  of  two,  however). 


3 . 0  Computation  of  The  Geomagnetic  Cutoff  Transmission  Function 

The  geomagnetic  cutoff  transmission  function  is  obtained  by  averaging 
the  geomagnetic  cutoff  around  the  spacecraft  orbit  so  as  to  obtain  a 
statement  of  the  probability  that  a  particle  can  reach  the  spacecraft  as 
a  function  of  Its  magnetic  rigidity.  The  resulting  transmission  function 
Is  used  to  modulate  the  particle  spectra  found  In  the  1 nterpl anetary 
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medium.  These  modulated  spectra  can,  In  turn,  be  used  to  estimate 
single-event  upset  rates. 

This  procedure  averages  over  the  large  variations  In  unset  rate  that 
are  to  be  expected  from  point  to  point  along  the  spacecraft's  orbit.  The 
only  justification  for  this  orbit-averaging  procedure  Is  that  the  mean 
time  between  upsets  Is  normally  found  to  be  long  compared  to  the 
spacecraft's  orbital  period.  For  each  spacecraft  design,  the  orbit 
average  upset  rate  should  be  checked  to  see  that  It  Is  no  more  than  a  one 
or  two  upsets  per  orbital  period.  If  the  predicted  rate  Is  higher,  then 
the  orbit  averaging  procedure  described  here  Is  not  valid  and  the 
variation  in  upset  rate  around  the  orbit  must  be  considered. 

By  using  the  vertical  cutoff  as  typical  of  all  directions,  employing 
eq.  (1)  for  altitude  extrapolation,  and  neglecting  the  earth's  shadow,  we 
have  simplified  the  geomagnetic  cutoff  transmission  calculation.  All 
that  is  necessary  now  to  obtain  the  orbit  -  averaged  geomagnetic 
transmission  to  the  skin  of  a  spacecraft  Is  to  compute  the  spacecraft's 
location  In  small  steps  around  the  orbit  and  Interpolate/extrapolate  the 
cutoff  value  to  each  point.  The  resulting  cutoffs  are  binned  Into  200 
bins  between  0  and  20  GeV/ec.  The  spacecraft  must  be  followed  for  a 
couple  of  days  so  that  Its  orbit  Is  allowed  to  precess  around  the  earth, 
sampling  the  worldwide  pattern  of  geomagnetic  cutoff.  Finally  the  binned 
cutoffs  are  normalized  to  obtain  the  transmission  function.  This 
procedure  Is  similar  to  the  one  reported  by  Heinrich  and  Spill  (1979). 


4.0  Computation  of  the  Spacecraft's  Position 

For  spacecraft  In  circular  orbits,  the  geographic  latitude  ,e,  of  the 
spacecraft  at  some  time,  t.  Is  given  by: 

o  «  cos-*  [sine0  s1n(ut  +  v0)],  (4) 

where  e0  Is  the  orbital  Inclination, 

t  Is  the  time, 

u  Is  the  satellite  angular  velocity  (radians/sec), 


%  Is  the  angle  between  the  spacecraft's  position  at  t  *  o  and 
the  point  on  the  equator  where  the  spacecraft  last  crossed  It  heading 
Into  the  northern  hemisphere.  This  angle  is  measured  about  the 
center  of  the  earth  (l.e.,  the  displacement  from  the  ascending  node). 

The  spacecraft  east  longitude,  I,  is: 

&  a  «f0  "  n  t  +  tan'l(Si)  +  tan"l(S2)  (5) 
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where 


cos  [1/2  (tr/2  -  b0)] 

cos  [1/2  ( ir/2  +  eQ)]  tan  [1/2  (n/2  -  ut  -  dQ)] 


,  =  sin  [1/2  (»/2  -  a  ) ]  , 

sin  [1/2U/2  +  e0)]  tan  [1/2  (t/2  -  «t  -  j*Q)] 

The  additional  symbols  in  these  equations  are  defined  as  follows: 
£1  is  the  angular  velocity  of  the  earth  (7.27X10-5  radians/sec) 


i0  is  the  longitude  of  the  point  on  the  equator  where  the 
spacecraft  last  crossed,  going  into  the  northern  hemisphere,  prior 
to  reaching  its  position  at  t  *  0  (i.e.  the  longitude  of  the 
ascending  node). 


For  a  circular  orbit,  the  spacecraft's  angular  velocity  is 
u  ■  1.24X10-3  [Re/(h+Re) ]3/2  radians/sec 


where 


Re  is  the  earth's  radius  (6371  km) 


h  is  the  spacecraft  altitude  in  km 

Now  we  need  to  compute  the  quantities  r,  x,  and  y  in  eq.  (1),  l.e. 
the  spacecraft's  position  and  the  particle's  arrival  direction  in  offset 
dipole  coordinates.  These  quantities  must  also  be  computed  for  the  four 
nearest  points  directly  below  the  spacecraft  on  the  world-wide  vertical 
cutoff  grid  at  20  km.  This  is  done  using  the  methods  of  Smart  and  Shea 
(1977).  With  these  quantities  determined  at  the  spacecraft's  position 
and  the  four  nearest  grid  points,  eq.  (1)  is  used  to  extrapolate  the 
cutoff  from  each  grid  point  to  the  spacecraft's  position.  These  four 
extrapolations  are  weighted  according  to  the  distance  between  the  point 
directly  beneath  the  spacecraft  on  the  20  km  altitude  world  grid  and 
each  grid  point,  then  they  are  averaged  to  obtain  the  estimate  of  the 
vertical  cutoff  at  the  spacecraft. 

Equations  (4)  through  (8)  were  used  to  compute  the  cutoff  at  the 
spacecraft's  position  at  200  points  per  orbit.  The  computations  were 
carried  out  for  two  days  of  orbits  to  obtain  good  convergence.  The 
computer  program  for  this  calculation  is  given  in  appendix  A  and  the 
geomagnetic  cutoff  tabulation  at  20  km  altitude  [taken  from  Shea  and 
Smart  (1975)]  is  given  in  appendix  B. 
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5.0  Comparison  with  Detailed  Calculations 


As  we  saw  in  section  2,  the  geomagnetic  cutoff  is  drastically 
different  at  different  points  around  a  spacecraft's  orbit.  It  is  also 
different  for  each  side  of  the  spacecraft  at  the  same  point  on  the 
orbit.  Finally  the  di rectlon-averaged  flux  never  reaches  the  value 
found  in  Interplanetary  space,  even  at  high  magnetic  rigidities.  This 
is  because  of  the  earth's  shadow  on  the  spacecraft. 

Drs.  D.  F.  Smart  and  M.  A.  Shea  of  Air  Force  Geophysical 
Laboratory  have  carried  out  detailed  calculations  of  the  geomagnetic 
cutoff  at  400  km  altitude  for  a  large  sampling  of  locations  and 
arrival  directions.  Recently,  we  have  used  these  results  to  make  a 
precise  calculation  of  the  geomagnetic  cutoff  transmission  function  to 
a  spacecraft  in  a  400  km  circular  orbit  at  several  inclinations  up  to 
50°.  These  results  take  into  account  the  cutoff  for  all  arrival 
directions  and  properly  account  for  the  earth's  shadow.  The  cutoff 
values  themselves  are  precise,  based  on  a  large  number  of  ray-tracing 
calculations  at  400  km.  Figures  2  and  3  show  the  transmission 
function,  averaged  over  all  arrival  directions,  at  two  different 
points  for  a  50°  inclination  orbit  at  this  altitude.  Figure  2  is  for 
10°  latitude  and  120°  longitude.  The  position  is  close  to  the 
highest  cutoff  to  be  found  around  any  orbit-  Figure  3  shows  the 
transmission  function  at  -50°  latitude  (southern  hemisphere)  and 
120°  longitude.  This  position  is  close  to  the  lowest  cutoff  to  be 
found  around  this  orbit.  The  cutoff  varies  from  -  0.4  GV  to  over 
10  GV  between  these  two  positions.  In  figure  3,  the  transmission 
rises  rapidly  to  its  limiting  value  (limited  by  the  earth's  shadow). 

By  contrast,  In  figure  2  the  transmission  does  not  reach  the  limiting 
value  until  nearly  100  GV.  This  shows  how  much  the  cutoff  depends  on 
the  particle's  arrival  direction  at  low  latitudes.  This  effect  can 
also  be  seen  in  eq.  (1). 

Figure  4a  shows  the  precise  geomagnetic  cutoff  transmission 
function,  averaged  over  all  arrival  directions  and  around  a  50° 
Inclination  circular  orbit  at  400  km  (curve  a).  If  this  function  Is 
multiplied  by  a  particle  spectrum  (in  particles/m2,  ster.  sec. 

MeV/u)  In  the  Interplanetary  space,  one  obtains  the  particle  spectrum 
at  the  skin  of  the  spacecraft  In  the  same  units. 

We  now  compare  this  precise  calculation  with  the  approximate 
method  described  in  sections  2.0,  3.0,  and  4.0.  Here  we  have 
neglected  the  change  In  cutoff  with  arrival  direction,  always  using 
the  vertical  cutoff.  We  have  also  employed  Stormer's  theory  to 
extrapolate  the  cutoff  from  precise  cutoff  values  at  20  km.  Figure  4 
shows  the  approximate  geomagnetic  cutoff  transmission  function 
neglecting  the  earth's  shadow  (curve  b)  and  the  same  function 
corrected  (approximately)  for  the  earth's  shadow  using  eq.  3  (curve  c). 

Figures  5,  6,  7,  8,  and  9  show  the  precise  transmission  function 
for  a  400  km  circular  orbit  at  40°,  30°,  20°,  10°  and  0° 
inclination  respectively.  Here  as  in  figures  1,  2,  and  3,  the 
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transmission  value  Is  expressed  as  a  solid  angle  (4*  correspondl ng  to 
a  transmission  factor  of  1.0). 


6 . 0  Computed  Geomagnetic  Cutoff  Transmission  Functions. 

The  geomagnetic  cutoff  transmission  function  has  been  computed,  by 
the  approximate  method,  for  circular  orbits  of  30°,  40°,  50°, 

60°,  70°,  80°,  and  90°  Inclination.  The  results  for  223  km 
altitude  are  given  In  figures  10.  Figure  11  gives  the  results  at  426 
km,  figure  12  at  741  km  and  figure  13  at  1111  km.  In  each  of  these 
figures,  the  lowest  cutoff  portion  of  the  transmission  function  for 
the  80°  and  90°  Inclination  orbits  Is  quite  uncertain.  This  Is 
because  the  geomagnetic  cutoff  near  the  earth's  poles  Is  unstable  and 
poorly  known.  Figure  14  compares  the  geomagnetic  cutoff  transmission 
function  to  be  expected  during  a  magnetic  storm,  with  the  same 
function,  taken  from  figure  13,  for  geomagnetical ly  quiet  times.  Both 
functions  are  for  a  60°  Inclination  orbit  at  1111  km.  The  earth's 
shadow  has  been  neglected  In  all  the  calculations  presented  In  this 
section. 


7.0  Use  of  the  Geomagnetic  Cutoff  Transmission  Function. 

To  obtain  the  differential  energy  spectra  for  the  various  nuclei 
reaching  the  skin  of  a  spacecraft  from  outside  the  magnetosphere,  just 
multiply  the  flux  In  the  Interplanetary  medium  (taken  from  Adams  et 
al.,  1981)  by  the  transmission  function.  To  do  this,  the  magnetic 
rigidity,  P  (In  GeV/ec)  must  be  computed  for  each  particle  energy,  E 
(In  MeV/u)  1 .e. : 

P  -  j  [(E/1000)2  +  1. 86X10“ 3E]1/2  (9) 

and  then  used  to  look  up  the  geomagnetic  cutoff  transmission.  Note  that  A 
Is  the  atomic  mass  and  Z  Is  the  particle  charge  (In  electron  charges). 

From  eq.  (9),  we  can  see  that  the  geomagnetic  transmission  at  any  particle 
energy  depends  on  the  particle's  charge.  If  an  Ion  Is  fully  stripped, 
then  A/Z  -  2,  however,  If  the  Ion  is  only  singly  Ionized,  A/Z  *  A.  Figure 
15  shows  the  oxygen  spectrum  In  the  Interplanetary  medium  (curve  a)  at  the 
skin  of  a  spacecraft  In  a  30°  Inclination  orbit  at  223  km  assuming  the 
oxygen  Is  fully  Ionized  (curve  b)  and  that  the  anomalous  component  (see 
Adams  et  al.,  1981)  of  oxygen  Is  singly-ionized  (curve  c). 


8.0  Conclusion 

This  report  presents  a  method  of  computing  the  geomagnetic  cutoff 
transmission  function  that  makes  it  possible  to  compute  the  orbit  - 
averaged  charged  particle  spectra  the  skin  of  any  spacecraft  in  a  circular 
orbl t. 


The  results  show  that  these  charged  particle  spectra  depend  critically 
on  the  particle's  charge  and  the  spacecraft's  orbital  inclination.  The 
spectra  also  depend,  to  a  lesser  degree  on  the  orbital  altitude. 

The  geomagnetic  cutoff  and  hence  the  charged  particle  spectra  change 
drastically  from  point  to  point  around  the  spacecraft  orbit.  At  each 
point  on  the  orbit,  the  particle  spectra  are  different  on  each  side  of  the 
spacecraft.  Also,  the  flux  striking  one  side  of  the  spacecraft  will  be 
reduced  by  the  cosmic  ray  shadow  of  the  earth.  All  of  these  points  have 
been  neglected  or  averaged  in  the  results  presented  here. 

The  points  discussed  above  have  two  practical  implications  for  upset 
rate  predictions:  1)  When  the  mean  time  between  upsets  Is  smaller  than 
the  orbital  period,  the  observed  upset  rate  at  the  low  cutoff  points  in 
the  orbit  can  be  dramatically  higher  than  predicted  by  this  procedure, 
especially  during  a  solar  flare.  2)  For  a  spacecraft  that  Is  stabilized, 
the  observed  upset  rate  will  be  dependent  on  the  physical  location  of  the 
microelectronic  components  within  the  spacecraft.  In  general,  the  lowest 
upset  rates  will  be  observed  when  the  electronics  are  located  on  the  side 
of  the  spacecraft  that  is  always  In  the  shadow  of  the  earth.  The  best 
location  will  differ  from  mission  to  mission  depending  on  the  details  of 
the  orbit  and  spacecraft  design.  If  the  trapped  radiation  is  the  dominant 
cause  of  upsets,  then  the  optimum  location  will  again  be  different. 
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EXPOSURE  FACTOR 


Figure  1.  The  precise  geomagnetic  cutoff  transmission  function, 
averaged  over  all  arrival  directions,  for  a  point  at  20°  latitude, 

0°  longitude  and  400  km  altitude.  The  arrow  shows  the  value  of  the 
the  vertical  geomagnetic  cutoff  at  this  location,  extrapolated  up  from 
20  km  using  eq.  (15.  This  function  Is  multiplied  by  the  particle 
spectra  In  Interplanetary  space  to  obtain  the  omnldi recti onal  spectra 
at  this  point.  The  exposure  factor  is  the  solid  angle  through  which 
cosmic  rays  may  arrive  at  this  point. 


RIGIDITY  (GV) 

Figure  2.  The  precise  geomagnetic  cutoff  transmission  function  at 
100  latitude,  120°  longitude  and  400  km  altitude.  This  Is  near 
the  location  of  the  highest  cutoff  experienced  by  any  spacecraft. 


EXPOSURE  FACTOR 


Figure  3.  The  precise  geomagnetic  cutoff  transmission  function  at 
-50°  latitude  (in  the  southern  hemisphere),  120°  longitude  and  400 
km  altitude.  This  is  near  the  location  of  the  lowest  cutoff 
experienced  by  a  spacecraft  in  a  50°  inclination  orbit. 


EXPOSURE  FACTOR 


Figure  4.  The  geomagnetic  cutoff  transmission  function  for  a 
spacecraft  in  a  50°  inclination  circular  orbit  at  400  km  altitude. 

The  cutoff  has  been  averaged  over  all  arrival  directions  and  over  the 
spacecraft's  orbit.  Curve  (a)  is  a  precise  calculation  based  on  a 
large  number  of  ray  tracing  calculations;  curve  (b)  is  an  approximate 
calculation  which  neglects  the  earth's  shadow;  and  curve  (c)  is  an 
approximate  calculation  which  takes  the  earth's  shadow  into  account  in 
an  approximate  way.  The  vertical  axis  here  runs  from  0.0  to  1.0, 
corresponding  to  the  range  0.0  to  4u,  in  the  previous  figures. 
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RIGIDITY  (GV) 

Figure  6.  The  same  as  figure  5,  but  for  30°  Inclination 


EXPOSURE  FACTOR 


0°  INCLINATION 
400  km  ALTITUDE 


1.0  10.0 

RIGIDITY  IGV) 


Figure  9.  The  same  as  figure  5,  but  for  QO  Inclination 


EXPOSURE  FACTOR 


RIGIDITY  (GV) 

Figure  10.  The  approximate  geomagnetic  cutoff  transmission  function 
for  a  circular  orbit  at  223  km  altitude.  The  earth's  shadow  has  been 
neglected.  The  curves  are  for  orbital  In'  '1  nations  of  90°,  80°, 

70°,  60°,  50°,  40°,  and  30°  as  shown.  To  obtain  the 
omnidirectional  spectra  at  the  skin  of  the  spacecraft  (in  particles 
/m2  sec.  MeV/u)  multiply  this  function  by  4ir  and  then  by  the  spectra 
(In  particles/m2,  ster.  sec.  MeV/u).  The  cutoff  probability, 
plotted  here  Is  the  exposure  factor  (In  figure  1)  divided  by  4ir. 
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EXPOSURE  FACTOR 


EXPOSURE  FACTOR 


Figure  14.  The  effect  of  a  magnetic  storm  on  the  geomagnetic  cutoff 
transmission  function  for  a  50°  inclination  orbit  at  1111  km:  The  lower 
curve  is  for  geomagnetical ly  quiet  times  and  Is  the  same  as  the  60° 
curve  In  figure  13.  The  upper  curve  is  for  geomagnetical ly  disturbed 
times  such  as  usually  exists  during  a  solar  flare  particle  event. 
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Figure  15.  The  oxygen  spectrum  at  solar  minimum.  Curve  (a)  is  the 
spectrum  as  seen  in  the  interplanetary  medium  near  earth.  Curve  (b)  is 
the  spectrum  at  the  skin  of  a  spacecraft  in  a  223  km  orbit  at  30° 
inclination  assuming  the  anomalous  component  is  fully-ionized.  Curve  (c) 
Is  the  same  as  curve  (b),  but  assuming  the  anomalous  component  is  singly 
Ionized. 


Appendix  A:  The  GEOMAG  Program 


The  following  program  Is  written  for  a  VAX  11/780  system.  A  similar 
version  has  been  used  on  a  PDP-10  computer.  The  program  was  written  to 
find  the  table  of  world-wide  vertical  cutoffs  In  a  file  called  CUTOFF.DAT 
(appendix  3).  The  geomagnetic  cutoff  transmission  function  is  tabulated 
and  written  on  an  output  file  called  GTRANS.DAT. 

This  program  Interacts  with  the  user  by  asking  a  series  of  questions 
to  obtain  the  data  needed  to  generate  the  GTRANS.DAT  file.  Following  the 
program  listing,  there  Is  an  example  of  a  run  and  a  listing  of  the 
resulting  GTRANS.DAT  file.  This  example  should  be  used  to  debug  the 
program  on  a  new  computer. 

The  program  contains  provisions  for  specifying  the  longitude  of  the 
ascending  mode  and  the  displacement  from  the  ascending  mode.  The  values 
set  by  the  program  are  zero  In  both  cases.  Provisions  are  also  made  for 
typing  out  the  cutoff  and  spacecraft  location  at  each  time  step  and  to 
compute  the  cutoff  In  any  direction  other  than  vertical.  The  function 
STORMER  contains  a  common  block  named  KARL  that  gives  access  to  the 
coordinates  of  the  spacecraft  and  the  particle's  arrival  In  offset  dipole 
coordinates. 
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6200 
6300 
6400 
6300 
6600 
6700 
6800 
6900 
7000 
7100 
7200 
7300 
7400 
7300 
7600 
7700 
7  BOO 
7900 
8000 
8100 
8200 
8300 
8400 
8300 
8600 
8700 
8800 
8900 
9000 
9100 
9200 
9300 
9400 
9300 
9600 
9700 
9800 
9900 
10000 
10100 
10200 
10300 
10400 
10300 
10600 
10700 
10BOO 
109u0 
11000 
11100 
11200 
11300 
11400 
11000 
11600 
11700 
11800 
11900 
12000 
12100 
12200 


C  READ  IN  THE  TABLE  OF  WORLD  WIDE  VERTICAL 

C  QEOMAQNETIC  CUTOFFS  AT  20KM  ALTITUDE 

C  THESE  DATA  ARE  TABULATED  EVERY  3  DE0REE8  IN  LATITUDE 

C  AND  EVERY  IS  DEOREEB  IN  LONGITUDE 

C  THEY  ARE  TAKEN  FROM  “M.  A.  SHEA  AND  D.  F.  SMART.  REPORT  NO. 

C  AFCRL-TR-73-0183.  HANSCOM  AFB. MASS.  1973" 

DO  20  1-1.  33 

II-34-I 

DO  20  J-l. 24 

READ <5,  403)  CUTOFF<II,J) 

C  THIS  ACCOUNTS  FOR  GEOMAGNETIC  CUTOFF  SUPRESSION  DURING 

C  LARGE  MAGNETIC  STORMS (FOLLOW I NO  ADAMS  ET  AL.  .  19B1) 

IF< ISTORM.  EQ,  1 )  CUTOFFt II, J ) -CUTOFF < 1 1 . J>*<1.  34*EXP( 

#  -CUTOFF ( 1 1,  J) /2.  9) ) 

20  CONTINUE 

C  AVERAGE  THE  CUTOFF  AROUND  THE  BO  DEGREE  LATITUDE  LINES, 

C  NORTH  AND  SOUTH 

CN-O. 

CS-O. 

DO  30  J-l. 24 
CN-CN+CUT0FFO3,  J) 

30  CS-C S+CUTOFF ( 1,  J) 

CN-CN/24. 

CS-CB/24. 

C  COMPUTE  THE  TOTAL  NUMBER  OF  STEPS  IN  TWO  DAYS  IF  WE  MAKE 

C  200  STEPS  PER  ORBIT 

JMAX-INT ( 400.  #W2/W1+1.  3> 

C  COMPUTE  THE  STEP  SIZE  IN  SECONDS 

STEP-2.  4PI/W2/200. 

C  COMPUTE  THE  VERTICAL  CUTOFF  AT  THE  SPACECRAFT 

r.  POSITION  FOR  EVERY  TIME  STEP 

DO  30  J-l.JMAX 
T1 -FLOAT (J-1)*STEP 
R 1-SIN <  THO ) #BIN ( W2*T 1+PBI ) 

THP— ACOS  <R 1 ) 

7LAT-90.  -180.  #THP/PI 

IF ( ZLAT.  QE.  80.  .  DR.  2LAT.  LE.  -80.  >  GO  TD  100 
C  COMPUTE  THE  SPACECRAFT  LONGITUDE 

RP-.  3* < PI/2.  +THO ) 

RM-.  3* < PI/2.  -THO) 

RF-.  3*<PI/2.  -W2#T1-PSI> 

IF  (  SIN<  RF ) .  EQ.  0.  >.  00  TO  63 

51- SIN(HM) *CQS(RF ) /SIN1RP  > /SIN(RF) 

52- C08<  RM) *CQS<RF ) /COS(RP ) /BIN(RF) 

SUM- ATAN  <  S 1 ) +ATAN ( 82 ) 

PHP— PHO-W 1 *T1 +SUM 
GO  TO  65 

63  PHP-P I +PHQ-W1 *T l 

63  IF  (PHP.  GE.  0.  .  AND.  PHP.  LT.  2.  #PI)  00  TO  70 

IF  (PHP.  LT.  0.  >  PHP-PHP+2.  *PI 
IFtPHP.  GE.  2.  #PI>  PHP-PHP-2.  *PI 
00  TO  63 
70  CONTINUE 

ZL0N-130,  *PHP/P I 

C  COMPUTE  THE  TABULAR  POSITION  OF  THE  VERTICAL  CUTOFF 

ZI-ZLAT/S.  +17. 

ZJ-ZLON/13.  +1. 

ILO-INT <  ZI ) 

IUP-ILD+1 

JLO-INT(ZU) 

JUP-JLO+1 
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12300 

12400 

12500 

12600 

12700 

12B00 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 

13000 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14000 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

15600 

15700 

15000 

15900 

16000 

16100 

16200 

16300 

16400 

16500 

16600 

16700 

16B00 

16900 

17000 

17100 

17200 

17300 

17400 

17500 

17600 

17700 

17000 

17900 

16000 

16100 

10200 

1B30O 


IF  (JUP.EQ.  25)  JUP-1 

C  INTERPOLATE  THE  VERTICAL  CUTOFF  TO  THE  EXACT  LOCATION 

C  OF  THE  SPACECRAFT  U8INQ  STDRMER  THEORY. 

DI-ZI-FLOATULO) 

DJ“Z J-FLOAT  < JLO) 

X0RB»(6371.  2+ALT 1/6371.  2 
RORD-  6391.  2/6371.  2 
ZLONLO" ( JLO- 1 1*15. 

ZLONUP-ZLONLQ-t-15. 

ZLATLO- IL0#5.  -B5. 

ZLATUP-ZLATLO+5. 

SC-STORMER<ZLAT,  ZLON, XORB,  AZ,  ZE) 

SCLL-STQRMER  <  ZLATLO. ZLONLO.  RORD,  AZO,  ZEO) 

SCUL-STQRMER <  ZLATUP,  ZLONLO.  RORD.  AZO.  ZEO) 

SCLU-STORMER  <  ZLATLO.  ZLONUP. RORD. AZO.  ZEO  > 

BCUU-STORMER ( ZLATUP.  ZLONUP.  RORD.  AZO.  ZEO) 

Yi»SC*CUTOFF< ILO.  ULOl/BCLL 
Y2-  SC*CUTOFF< IUP,  JLOl/SCUL 
Y3-  SC*CUTQFF  < ILO. JUP ) /SCLU 
Y4-  SC*CUTOFF< IUP. UUP ) /SCUU 

C»<1. -D1 ) #< 1.  -DJ)#Y1+<1.  -DI >*DJ*Y3+DI*< 1.  -DJ)*Y2+DI*DJ*Y4 
00  TO  200 

100  IFIZLAT,  LE.  -SO.  >00  TO  110 

SC-STORMER<90.  >  0. . XORB, AZ. ZE> 

SC  X-STORMER  ( 90.  .  0.  ,  RORD,  AZQ,  ZEO) 

C-BC*CN/SCX 
00  TO  200 

110  SC-STORMER < -90,  ,0.  ,  XORB.AZ.  ZE) 

BCX«ST0RMER(-90.  ,0.  , RORD, AZO. ZEO > 

C-SC#CS/SCX 

C  HISTOORAM  THE  CUTOFFB 

200  IDEX-1NT(C#10,  )  +  l 

MAT ( IDEX )-MAT  < IDEX )  +  1 

C  #*##***«**##*•*  **•**#******##**#**##*****###**#*#***##*■*#**•*#*■*•**#*# 

C  TO  PRINT  OUT  THE  LONOITUDE,  LATITUDE,  AND  CUTOFF  AT  EACH 

C  POINT  ARROUND  THE  SPACECRAFT'S  ORBIT.  JUST  ENABLE  THE  FOLLOWING 

C  STATEMENTS  BY  REMOVINO  THE  COMMENT  C'S. 

C##*»*»*«**»*»»*******#***»**»#*«****#**#****ft*#**»******* ***#*##*#*#* 

C  TYPE  435,  Zl  ON,  ZLAT,  C 

50  CONTINUE 

CMAT-O. 

C  SAVE  THE  THRESHOLD 

JSAV-0 

DO  300  J- 1,200 

C  CONVER  THE  HISTOORAM  TO  TRANSMISSION 

CHAT  “FLOAT  <  MAT  <  J ) )  /FLOAT  ( JMAX  >  t-CMAT 
IF <  (JSAV  EQ.  0).  AND.  <CMAT.  OT.  0)  >  JSAV-J+1 
CF ( J+l ) "FLOAT  <  J ) / 10. 

300  T<  J+U-CMAT 
CFC l)-0. 0 
T<l>-0.  0 

C  ADJUST  THE  HIGHEST  ZERO  TRANSMISSION  VALUE  TO  BE 

C  JUST  AT  THE  CUTOFF  THRESHOLD 

CF  <  JSAV- 1  > -CF  ( JSAV )  -T (  JSAV  >  # ( CF  <  JSAV+i  >  -CF  <  JSAV )  >  / 

1  <  T  <  J8 AV+ 1 ) -T  <  JSAV ) ) 

C  THIS  IS  A  CORRECTION  FOR  THE  EARTH'S  SHADOW  ON  THE  SPACECRAFT 

C  ACCORDINO  TO  SIMPLE  GEOMETRICAL  OPTICS. 

IF < ISHADOW,  NE. 1 )  00  TO  390 
DO  380  J- 1,200 

T< J)-T<J>*<1. -0.  5* ( 1 ,  -< (6371.  2+ALT>»#2.  -<6371.  2)*#2)**.  5/ 

*  <6371.  2+ALT) ) ) 
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18400  380  CONTINUE 

18500  390  00  400  J- 1.200 

18600  400  WRITE!6.  410JCF!u> , T! J> 

18700  403  FORMAT ( 17X. F6. 2) 

18800  410  FORMAT <  3X.  F6.  3.  5X< FB,  6) 

18900  411  FORMAT < IX/  'DO  YOU  WANT  TO  INCLUDE  THE  EFFECT  OF  THE  SHADOW'./, 

19000  *  '  OF  THE  EARTH?  <1  FOR  YESi  0  FOR  NO)  '.*> 

19100  412  FORMAT! IX,  'ENTER  THE  MAONETIC  WEATHER  CONTITION:  1  FOR  STORMYi 

19200  *  0  FOR  QUIET  '. 4) 

19300  413  FORMAT! IX.  'ENTER  ORBITAL  INCLINATION  IN  DECREES  '*> 

19400  420  FORMAT (IX.  'ENTER  LONGITUDE  OF  A8CENDXNQ  NODE  IN  DEGREES  '«> 

19300  423  FORMAT! IX.  'ENTER  DISPLACEMENT  FROM  NODE  IN  DEGREES  '*) 

19600  430  FORMAT! IX.  'ENTER  ALTITUDE  IN  KILOMETERS  '*> 

19700  433  FORMAT!'  LONOITUDE  «  '.F6.  2.  '.  LATITUDE  -  '.F3.  2.  '  CUTOFF  -  '. 

19000  #  F4.  3) 

19900  END 
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FUNCTION  STORMER ( OCLATD, OCLOkD, ROC. AZ. ZE> 

COMMON /KARL /RED.  EDLAT .  AZM.  ZEM.  GAMMA 

THIS  FUNCTION  TRANSFORMS  A  OEOQRAPHIC  LOCATION  AND  ARRIVAL 
DIRECTION  INTO  OFFBET  DIPOLE  COORDINATES.  THEN  COMPUTES  THE 
STORMER  CUTOFF  IN  OV  AND  RETURNS  THE  RESULT.  THE  OFFSET  DIPOLE 
COORDINATES  ARE  AVAILABLE  IN  THE  COMMON  BLOCK  /KARL/. 


GCLATD  IS  GEOCENTRIC  LATITUDE  IN  DEGREES 

OCLQNQ  IS  GEOCENTRIC  LONGITUDE  IN  DEGREES 

ROC  IS  RADIAL  DISTANCE  FROM  QEOCENTER  IN  EARTH  RADI  I 

AZ  IS  GEOGRAPHIC  AZIMUTH 

ZE  IS  GEOGRAPHIC  ZENITH 

RED  IS  RADIAL  DISTANCE  FROM  OFFSET  DIPOLE  POSITION  IN 

EARTH  RADII 

EDLAT  IS  THE  GEOMAONETIC  LATITUDE  IN  OFFSET  DIPOLE  COORDINATES 
AZM  IS  GEOMAGNETIC  AZIMUTH  IN  OFFSET  DIPOLE  COORDINATES 

ZEM  IS  GEOMAGNETIC  ZENITH  IN  OFFSET  DIPOLE  COORDINATES 

GAMMA  IB  GAMMA  ANGLE  MEASURED  FROM  MAGNETIC  EAST 

IF<  JDAT A. EQ.  77 )  00  TO  10 

PI  -  ACOB < - 1 . 0 ) 

RAD  -  ISO. O/PI 
PI 02  -  PI/2. 0 
TUOPI  ■  PI*2.  0 

SATA  ERAD,  THETAD.  PHID.  R1KM.  TH1DEG.  PH1DE0/6371. 2, 

1  11.  4354, -290.  2392,  450  2586, 72.  8278,  148.  7753/ 

NOPT  ■  0 

SORTS  -  SORT (3.0) 

ENTER  GEOMAGNETIC  DATA.  IORF  1975 


DT  -  ! 

SEE  JOR, 

5.  0 

Bl. 

5163,  1976 

DT  IB  NUMBER  OF  YEARS  SINCE  1975 

001  - 

-301S6.  0 

4- 

25.  6*DT 

002  - 

-1898.  0 

- 

24.  9*DT 

on  - 

-2036.  0 

♦ 

10.  0*DT 

012  - 

2997.  0 

+ 

0.  7*DT 

HI  1  - 

5735.  0 

- 

10.  2*DT 

H12  - 

-2124,  0 

- 

3.  0#DT 

022  - 

1551.  0 

4.  3*DT 

H22  - 

-37.  0 

- 

IB.  9#DT 

IF  (NOPT.  EO.  1 )  PRINT  1000,  001.  002,  Oil.  012.  022,  HU.  H12,  H22 
C  COMPUTE  POSITION  OF  OFFSET  DIPOLE 

HO  -  SQRT(  001*001+01 1*01 1+H1 1  *W1 1 > 

HOBO  ■  H0*H0 

ELO  -  2  0*001  *002+ (Q'l  1*G12+Hl  1*H12>  *SQRT3 
ELI  ■  ~Qll*Q02+<G01*Q12+011*022+Htl*H22>*SQRT3 
EL2  -  -HI 1*002+ <  Q01*H12-H1 1*022+01 1*H22)*SQRT3 
E  -  <  EL0*G01+EL1 *01 1+EL2*H1 1 ) *4,  0*H0SQ 
E  ■  <  EL 0*00 1+EL 1*01 1+EL2*H1 1 ) / <  4.  0*H0BQ) 

IF<NOPT. EQ.  1)  PRINT  1011.  ELO,  ELI,  EL2,  E.  HO 
1011  FORMAT ( 1H  ,  SE15,  5) 

XEDFOC  »  ERAD*<ELl-011*E>/<3,  0*H0S0> 

XEDFOC  *  <ELl-011*E)/(3.  0*H08Q) 

YEDFOC  *  ERAD* ( EL2-H1 1 *E ) / ( 3.  0*H0SG ) 

YEDFOC  -  < EL2-H1 1 *E ) / ( 3,  0*H0BQ) 

ZEDFOC  •  ERAD* (EL0-001 *E > / ( 3.  0*H0SQ) 

ZEDFOC  •  < EL 0-001 *C ) / t  3.  0*H0SQ> 

REDFOC  *  SORT <  XEDFOC*XEDFOC+YEDFOC*YEDFOC+ZEDFOC*ZtDFOC ) 

IF(NOPT.  EQ.  UPRINT  3001,  XEDFOC,  YEDFOC.  ZEDFOC,  REDFOC 
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3001  FORMAT  < 1H  .  4F10. 4.  3X,  'XEDFQC,  YEDFGC,  ZEDFQC,  REDFOC ' ) 

1000  FORMAT  ( 1H  .  10F13.  5) 

1010  F0RMAT(1H0.  9F19,  9/1H  , BF19.  9) 

THETA  -  THETAD/RAD 
PHI  -  PHID/RAD 
CP  -  C03(PHI > 

SP  -  SIN(PHI) 

ST  -  SIN(THETA) 

CT  -  CQB(THETA) 

CPCT  -  CP*CT 
CPST  -  CP#ST 
SPOT  -  SP*CT 
SPST  -  8P*ST 
R1ER  -  H1HM/ERAD 
TH1RAD  -  TH1DE0/RAD 
PH1RAD  -  PH1DEQ/RAD 

IF < NOPT , EQ.  DPRINT  1000.  RIKM,  TH1DE0.  PH1DEQ,  R1ER.  TH1RAD,  PH1RAD 

XQMED  ■  XEDFQC *CPCT  -YEDFQC»BPCT  -ZEDF«C*BT 

YOMED  -  XEDFQC »BP  ♦YEDFQC^CP 

ZONED  «•  XEDFQCaCPBT  -VEDFQC*BPBT  +ZEDFQC»CT 

IF <NOPT. EO.  1 ) PRINT  3002.  XQMED.  YQMED,  ZQMED 

3002  FORMAT < 1H  .  3F10  4.  13X.  'XQMED.  YQMED,  ZQMED') 

IF (NOPT. EQ.  DPRINT  1010.  CP.  BP,  CT,  BT.  CPCT.  CPST,  SPCT .  8PST 
JDATA  -  77 

10  CONTINUE 

ITERATE  TO  FIND  COORDINATES  OF  OFFBET  NORTH  DIPOLE  AT  ANY 
LATITUDE 

FIRBT  QUES3  FIND  OFFSET  NORTH  DIPOLE  AT  DISTANCE  ROC 
ZDEDNP  -  ROC 

100  XODNP  -  XQMED#CPCT  +  YQMEOmSP  +  ZDEDNP *CP6T 
YODNP  -  -XQMED*SPCT  +  YGMED*CP  -  ZDEDNP*SPBT 
ZODNP  -  -XOMED*BP  +  ZDEDNP#CT 

DODNP  -  SORT ( XODNP#XODNP  YODNP*YODNP  +  ZODNP*ZODNP> 

DIFLA  -  DODNP  -  RQC 

IF (ABB<  DIFLA)  -  l.OE-9)  ISO,  120.  110 
110  PHINOF  *  ATAN2 (YODNP,  XODNP )#R AD 

IF(PHINQF.  LT.  0.  0)  PHINOF  -  PHINOF  *  340. 0 
TNOF  -  -ACOS ( ZODNP /DODNP ) *RAD  *  90.0 

IF (NOPT. EQ.  DPRINT  4001 ,  ZDEDNP,  XODNP.  YODNP,  ZODNP, DODNP.  DIFLA. 

1  PHINOF,  TNOF 

4001  FORMAT  ( 1H  .  SX.  'ODC  0,  0,  'F7.  9,  '  -  QC  X,  Y,  Z  0F'3FB.  9, 

1  '  DODNP  -'F9.5.  '  DIF  OF  'F9.  4,  '  AT  LOND  LAT'FIO.4.  FB.4) 
ZDEDNP  -  ZDEDNP  -  DIFLA 
QO  TO  100 
120  CONTINUE 

PHINOF  -  AT AN2( YODNP, XODNP )«R AD 
IF<PHINOF.  LT.  0.  0)  PHINOF  -  PHINOF  +  340.0 
TNOF  -  -ACOS ( ZODNP/ DODNP )*R AD  +  90.0 

IF (NOPT. EQ.  DPRINT  4001,  ZDEDNP,  XODNP,  YODNP, ZODNP, DODNP, DIFLA, 

1  PHINOF,  TNOF 

BQCLATD  -  BIN(GCLATD/RAD) 

CQCLATD  -  CQB(QCLATD/RAD) 

8QCL0ND  -  8IN< QCLOND/RAD) 

CQCLOND  -  C08<  QCLOND/RAD ) 

C  QET  QEOCENTRIC  X  Y  Z  COORDINATES 

XQC  -  RQC*CQCLATD*CQCLOND 
YQC  -  RQC#CGCLATD*SQCLOND 
ZQC  «  RQC*SQCLATD 
QCT  -  (90.0  -  QCLATD)  /RAD 
SQCT  -  SIN(QCT) 

CQCT  -  COS (QCT) 
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FIND  X  Y  Z  IN  LOCAL  COORDINATES  OF  X-0,  Y-0,  Z 

THE  LOCAL  COORDINATE  Z  AXIS  PAS8ES  THRU  P 
THE  LOCAL  COORDINATE  X. Z  PLANE  CONTAINS  P 
OCROT  -  ATAN2  < YCC,  XQC ) 

IF<NOPT, EQ.  DPRINT  2001.  XQC.  YQC.  ZQC.  OCROT 

FORMAT < 1H  .  4F10,  4.  3X,  'XQC,  YOC.  7GC,  OCROT') 

BCCROT  -  SIN< OCROT) 

COCROT  -  COB ( OCROT ) 

XRL  -  XGC*CGCROT#CQCT  +  YOC»SOCROT*COCT  -  ZGC*SQCT 

YRL  -  -XGC*BGCROT  +  YCC*CQCROT 

ZRL  -  XGC*COCROT *SGCT  +  YQC*SGCROT*BQCT  ZGC#CGCT 


-  ZEDFGC*SGCT 


FORMATUH,  3F10.  4,  13X,  'XRL.  YRL.  ZRL') 

IF(NOPT. EO.  lJPRINT  2002.  XRL.  YRL.  ZRL 

DETERMINE  LOCATION  OF  OFFSET  DIPOLE  CENTER  IN  THEBE  SAME 

ROTATED  LOCAL  COORDINATES 

XEDRL  -  XEDFGC*CGCROT*CGCT  4-  YEDFOC*BGCRQT*CGCT  -  ZEDFQC*SGCT 
YEDRL  -  -XEDFQC#BOCROT  ♦  YEDFOC*COCROT 

ZEDRL  -  XEDFOC*CGCROT*BGCT  +  YEDFGC#SGCROT*SGCT  +  ZEDFOC*COCT 
IFINOPT.EQ.  DPRINT  3003.  XEDRL.  YEDRL.  ZEDRL 
FORMAT  <1H  .  3F10.  4.  13X,  'XEDRLM  YEDRL.  ZEDRL') 

XEDRTL  -  XEDRL 
YEDRTL  -  YEDRL 
ZEDRTL  -  ZEDRL  -  ZRL 

IF < NOPT.  EO.  DPRINT  2303.  XEDRTL.  YEDRTL.  ZEDRTL 


FORMAT  (1H  ,  3F10.  4, 


'XEDRTL.  YEDRTK.  ZEDRTL') 


TRANSLATE  TO  LOCAL  COORDINATE  SYSTEM  WITH  ORIGIN  AT  BUR FACE 
XRTL  -  XRL 
YRTL  -  YRL 
ZRTL  —ZRL 


YEDP  -  YRTL  +  YEDRL 
ZEDP  -  ZRTL  +  ZEDRL 

RED  -  SORT  <  XEDP*XEDP+YEDP*YEDP+ZEDP*ZEDP  > 

FORMAT  ( 1H  .  3F10. 4.  13X.  'XRTL.  YRTL,  ZRTL') 

IF < NOP T.  EG.  DPRINT  2302.  XRTL,  YRTL,  ZRTL 

EARTHS  SURFACE  AT  A  SPECIFIED  ALTITUDE 
POSITION  OF  OFFSET  NORTH  DIPOLE  IN  LOCAL  COORDINATE  8YBTEM 
XODNPR-  XQDNP*CGCROT*COCT  ♦  YQDNP-BQCRQT *CQC T  -  ZODNP*SQCT 

YODNPR-  -XODNP*SGCROT  +  YDDNP#CQCRQT 

ZODNPR-  X0DNP#CGCRQT#8QCT  +  YQDNP*90CR0T#BGCT  +  ZODNP#CQCT 

XQDNPT  -  XODNPR 
YODNPT  -  YODNPR 
ZODNPT  -  ZODNPR  -  ZRL 


ZODNP#CQCT 


XODNPR-  XODNP#CGCROT*COCT  ♦  YQDNP-BQCRQT *CQC T  -  ZODNP*SQCT 

YODNPR-  -XODNP*SGCROT  +  YODNP#CQCRQT 

ZODNPR-  X0DNP*CQCR0T#80CT  +  YQDNP#S0CR0T#80CT  +  ZODNP#CQCT 

XQDNPT  -  XODNPR 
YODNPT  -  YODNPR 
ZODNPT  -  ZODNPR  -  ZRL 

IF(NOPT.  EQ.  DPRINT  D03,  XODNPT,  YODNPT,  ZODNPT 
1  FORMAT ( IH  .  iOX,  'XODNPT  «'F10.  4, 2X.  'YODNPT  -'FtO. 4.2X,  'ZODNPT  -' 

1  F10. 4. 3X.  'OFFSET  N  DIPOLE  IN  LOCAL  COORDINATES') 

ROTM  -  ATAN2< YODNPT,  XQDNPT)  +  PI 

FIND  ANQLE  FROM  QEOQRAPHIC  NORTH 
NEGATIVE  -  ROTATION  FROM  QEOQRAPHIC  NP  CLOCKWISE 
POSITIVE  -  ROTATION  FROM  QEOQRAPHIC  NP  CCW 
SROTM  «  SIN(ROTM) 

CROTM  -  COS < ROTM) 

ROTMD  -  ROTM*RAD 

T  FORMAT  < IH  .  FI 5. D.  3X.  '  ROTM  IN  DEGREES  MEASURED  CCW  SO  -X 

l  WILL  POINT  TOWARD  OFFSET  NORTH  DIPOLE  AXIS') 

IFINOPT.  EQ.  DPRINT  2327.  ROTMD 

FTND  COMPONENTS  OF  UNIT  VECTOR  AT  ARBITARY  AZIMUTH  AND  ZENITH 
PLAZ  -  -AZ/RAD  +  PI 
TLZE  -  ZE/RAD 
SPLAZ  -  SIN(PLAZ) 

CPLAZ  -  COS (PLAZ) 


,  ”"Vf  ITT? 


r*-  ~  0 


o  o  o  o  o  o  o 


2003 


3001 

c 


3002 

c 

c 


1101 

1 

c 


c 

c 


c 

c 


STLZE  ■  SIN(TLZE) 

CTLZE  -  COS(TLZE) 

XLD  -  STLZE#CPLAZ 
YLD  -  STLZE#SPLAZ 
ZLD  -  CTLZE 

IF (NOPT. EQ.  D PRINT  2003,  XLD,  YLD,  ZLD,  AZ,  ZE 

FORMAT  <1H  ,  3F10.  4,  3X,  'UNIT  VECTOR  COMPOENTS  AT  AZ  it  ZE') 

FIND  COMPONENTS  OF  UNIT  VECTOR  IN  DIPOLE  RADIAL  COORDINATES 

ROTATE  AROUND  Y  AXIS  SO  -Z  AXIS  PASSES  THROUOH  XED.  0,  ZED 
NEW  VECTOR  IS  VA  -  ZRTL  «■  ZEDRTK  +  XEDRTL 
ANGLE  BETWEEN  VECTOR  FROM  POINT  LOCAL  ORIGIN  TO  GEOCENTER 
AND  VECTOR  FROM  POINT  LOCAL  ORIGIN  TO  XED,  0,  ZED 
CA  -  ZRTL*ZEDRTL/ ( ABS! ZRTL)*SORT<  ZEDRTL*ZEDRTL  +  XEDRTL*XEDRTL) > 

A  -  ACOS(CA) 

IF ( XEDRTL.  OT.  0.  0)  A  -  -A 
BA  -  BIN(A) 

ADEO  -  A*RAD 

lFINOPT.EQ.  DPRINT  100C,  CA.  A,  SA,  ADEG 

XLP  -  XI.D*CA  +  ZLD*9A 

YLP  -  YLD 

ZLP  •- XLD*SA  +  ZLD*CA 

IF<NOPT.  EO,  DPRINT  3001,  XLP,  YLP,  ZLP 

FORMAT! 1H  ,  3F104,  13X,  'XLP,  YLP,  ZLP  ') 

ROTATE  AROUND  X  PRIME  AXIS  SO  -Z  PABSES  THROUGH  XED,  YED,  ZED 
CB  -  ZRTL#ZEDRTLX  ( ABB<  ZRTL)*SQRT <  ZEDRTL#ZEDRTL  *  YEDRTL*YEDRTL ) ) 

B  «  ACOS!CB  > 

IF ( YEDRTL.  OT.  0.  0)  B  -  -B 

SB  ■  SIN(B) 

BDEQ  -  8»RAD' 

IF  (NOPT.  EG.  DPRINT  1000,  CB,  B,  SB,  BDEC 
XLPP  a  XLP 

YLPP  -  YLP#CB  +  2LP#SB 
ZLPP  — YLP*BB  ♦  ZLP#CB 

lFINOPT.EQ.  DPRINT  3002,  XLPP,  YLPP.  ZLPP 
FORMAT! 1H  ,  3F10. 4,  13X.  'XLPP,  YLPP.  ZLPP  ') 

ROTATE  AROUND  ZPP  AXIS  BO  -X  AXIS  PASSES  THROUOH  NORTH 

OFFSET  DIPOLE  AX  IB 

ZLDM  -  ZLPP 

XLDM  -  XLPP*CROTM  +  YLPP»8R0TM 
XLDM  -  XLPP*CROTM  -YLPP*SROTM 
YLDM  — XLPP#SROTM  +  YLPP#CROTM 
YLDM  -  XLPP*SROTM  +  Yl..PP*CROTM 

IF(NOPT.  EQ.  DPRINT  1101.  XLD,  YLD,  ZLD,  XLDM,  YLDM.  ZLDM 
FORMAT  < 1H  .  'UNIT  VECTOR  IN  LOCAL  COORDINATES  3F8.  3, 3X, 

'UNIT  VECTOR  IN  LOCAL  MAGNETIC  COORDINATES',  3F10.  3) 

FIND  AZUMITH  ANGLE  OF  UNIT  VECTOR  IN  LOCAL  DIPOLAR  RADIAL  COOR 
PAZM  -  ATAN2! YLDM,  XLDM) 

AZM  «  (PI  -  PAZM) *RAD 
IF ! AZM.  OT.  360.  0)  AZM  -  AZM  -  360.0 
ZEM  -  AC QB ( ZLDM ) #R  AD 
FIND  GAMMA  ANGLE 
GAMMA  ■  ACOB!  YLnt1)»RAD 
TRANSFORM  TO  OFFSET  DIPOLE  COORDINATES 
XED1-XQC-XEDFCC 
YED1-YQC-YEDFCC 
ZEDl-ZGC-ZEDFQC 

FIND  THE  Z  COORDINATE  IN  OFFBET  DIPOLE  COORDINATES 
ZED2">XED1#CPST-YED1*BPBT+ZED1  #CT 
FIND  THE  GEOMAGNETIC  LATITUDE 


EDLAT-R AD# ( P I Q2-ACQS ( ZEDS/RED) ) 

COSLDA-CQS < EDLAT/RAD ) 

STQRMER-AO.  #C0SI_DA##4.  !  <RED#RED* <  1,  4.SORTU.  -C0BLDA#*3.  *YLDM))*#a) 
RETURN 
END 


UsiVnam* I  ADAMS 
Paaaword  t 

Utlcoma  to  YAX/VMS  varaiori  V3.0 

******************************************************************** 
WELCOME  TO  THE  NRL  SPACE  SCIENCE  DIVISION  VAX  11/780  SYSTEM 
FOP  INFORMATION  ON  SYSTEM  ASSISTANCE!  t  TYPE  BYStMESSAGESASSIST.TXT 
************** **************** ************************************** 

*  R  GEOMAO 

DO  YOU  WANT  TO  INCLUDE  THE  EFFECT  OF  THE  SHADOW 
OF  THE  EARTH?  <1  FOR  YES>  0  FOR  NO)  0 

ENTER  THE  MAGNETIC  WEATHER  CONTITION!  1  FOR  STORMY  I  0  FOR  QUIET  0 
ENTER  ORBITAL  INCLINATION  IN  DEGREES  AO 
ENTER  ALTITUDE  IN  KILOMETERS  1111 

*  LO 

ADAMS  la**ad  out  at  29-N0V-19B2  1112ASS9.48 


RIGIDITY  TRANSMISSION 


0.  098  0.  000000 

0,  100  0.  00969B 

0.  800  0.  032824 

0.  300  0.  092406 

o'.  400  0.  069191 

0.  900  0.  084856 

0. 600  O.  09B697 

0.  700  0.  1 1 1899 

0,  800  O.  129699 

0  900  0.  138794 

1.000  0.191996 

1  100  0.  166396 

1.  200  0.  180903 

1.300  0.197901 

1.400  0.212794 

1.900  0.227194 

1.600  0.244125 

1.700  0.298299 

1 . 800  0.  270794 

1.900  0.  2B19B4 

2  000  0.293174 

2.100  0.309110 

2  200  0.319994 

2.  300  0.  324692 

2.  400  0.  336628 

2.  900  0.  344B34 

2.  600  0.  355278 

2.  700  0.  364603 

2l  BOO  0.  373555 

2.  900  O.  383626 

3.  000  0.  394256 

3.  100  0.  403767 

3.  200  0.  412346 

3  300  0.419806 

3.  400  0.  427B29 

3,  500  0.  435658 

3.600  0.4431  IB 

3.700  0,450019 

3.800  0.457106 

3.900  0.464192 


4.  000 

0.  469601 

4.  100 

0.  476315 

4.  200 

0.  481723 

4.  300 

0.  486759 

4.  400 

0.  492940 

4.  900 

0.  498135 

4.  600 

0.  503170 

4.  700 

0.  507460 

4.  BOO 

0.  5 12682 

4.  900 

0.  518277 

5.  000 

0.  524053 

9.  100 

0.  527975 

9.  200 

0.  531705 

5,  300 

0,  535B07 

5.  400 

0.  541775 

5.  900 

0.  547370 

5.  600 

0.  551660 

9,  700 

0  555949 

9.  BOO 

559866 

9.  900 

0.  564715 

6.  000 

0.  569191 

6.  100 

0.  973667 

6.  200 

0.  977769 

6.  300 

0.  9B2249 

6.  400 

0.  9B8213 

6.  500 

0.  992B76 

6.  600 

0.  999B60 

6.  700 

0.  999776 

6.  800 

0.  609371 

6.  900 

0.  610407 

7,  000 

0.  619069 

7.  100 

0.  61991B 

7.  200 

0.  624  900 

7.  300 

0.  629096 

7.  400 

0.  634B3B 

7.  500 

0.  63B199 

7.  600 

0.  644936 

7.  700 

0.  64B266 

7.  800 

0.  694900 

7.  900 

0.  660019 

B.  000 

0.  663745 

B.  100 

0.  670459 

B.  200 

0.  675121 

8.  300 

0.  679411 

8.  400 

0.  6B9379 

8.  900 

0.  692839 

8.  600 

0.  697 12B 

8.  700 

0.  702909 

8.  800 

0.  709250 

8.  900 

0.  713726 

9.  000 

0.  723051 

9.  100 

0.  729952 

9.  200 

0.  739733 

9.  300 

0.  743966 

9.  400 

0.  791212 

9.  900 

0.  761696 

9.  600 

0.  771941 

9.  700 

0.  7B3103 

9.  800 

0.  799789 

9.  900 

0.  806416 

10.  000 

0.  816673 

10.  100 

0.  829065 

10. 200 

0.  834950 

10.  300 

0.  845207 

10.  400 

0.  855091 

10. 900 

0.  862365 

10, 600 

0.  875233 

10. 700 

0.  BB41B5 

10. BOO 

0.  B9201B 

10. 900 

0.  900970 

11. 000 

0.  910481 

11.  100 

0.  917941 

1 1 ,  200 

Q.  926706 

11.  300 

0.  934166 

11.  400 

0.  9393BB 

11. 900 

0.  946475 

11.  600 

0,  950951 

11.  700 

0.  956919 

11.  BOO 

0.  960B35 

1 1 .  900 

0.  966057 

12.  000 

0.  970906 

12.  100 

0.  974263 

34 


"Vgfanx 


is.  aoo 

0  977993 

18. 300 

1.  OOOOOO 

ia.  300 

0.  983215 

18. 400 

1,  OOOOOO 

IS.  400 

0.  986372 

18.  SOO 

1.  OOOOOO 

12.  900 

0.  989742 

18. 600 

1.  OOOOOO 

12.  600 

0.  9936S9 

18.  700 

1.  oooooo 

12.  700 

0.  996496 

is.  aoo 

1.  oooooo 

12,  800 

1.  OOOOOO 

18.  900 

1.  OOOOOO 

12,  900 

1.  OOOOOO 

19. 000 

1.  oooooo 

13,  000 

1,  OOOOOO 

19. 100 

l.  oooooo 

13.  100 

t.  OOOOOO 

19.  200 

1.  oooooo 

13.  200 

1.  OOOOOO 

19.  300 

1.  oooooo 

13.  300 

1.  OOOOOO 

19.  400 

1.  oooooo 

13,  400 

1.  OOOOOO 

19,  SOO 

1.  oooooo 

13.  SOO 

1.  OOOOOO 

19.  600 

1.  oooooo 

13.  600 

1,  OOOOOO 

19.  700 

1.  oooooo 

13.  700 

l.  OOOOOO 

19.  BOO 

1.  oooooo 

13.  800 

1.  OOOOOO 

19. 900 

1.  oooooo 

13.  900 

1.  OOOOOO 

14.  000 

1.  oooooo 

14.  100 

1.  oooooo 

14.  200 

1.  oooooo 

14. 300 

1.  oooooo 

14.  400 

1.  oooooo 

14.  900 

1.  oooooo 

14.  600 

1.  oooooo 

14.  700 

1.  oooooo 

14.  800 

1.  oooooo 

14. 900 

1.  oooooo 

IS,  000 

1.  oooooo 

IS,  too 

1.  oooooo 

IS.  200 

1.  oooooo 

IS.  300 

1.  oooooo 

IS.  400  1.  000000 

IS.  900  1.  000000 

IS.  600  1.  OOOOOO 

IS.  700  t.  OOOOOO 

is.  aoo  l.oooooo 

15.  900  t.  OOOOOO 

16.  000  1.  OOOOOO 

16.  100  1.  OOOOOO 

16.  aoo  l.  oooooo 

16.  300  1.  OOOOOO 

16.  400  1.  OOOOOO 

16.  SOO  1.  oooooo 

16.600  1.000000 

16.  700  1.  OOOOOO 

16.  800  1.  OOOOOO 

16.  900  1.  OOOOOO 

17.000  1.000000 

17.  100  1.  OOOOOO 

17.  aoo  i.  oooooo 

17.  300  1.  OOOOOO 

17.  400  1.  OOOOOO 

17.  SOO  1.  OOOOOO 

17,  600  1.  OOOOOO 

17.  700  1.  OOOOOO 

17.  800  1.  OOOOOO 

17.  900  1.  OOOOOO 

18.  000  1.  OOOOOO 

18.  100  l.  OOOOOO 

ia.  200  t.  oooooo 


Appendix  B!  Table  of  Vertical  Cutoffs 


The  following  data  table  Is  Input  to  the  GEOMAG  program  as 
CUTOFF.DAT.  It  consists  of  the  vertical  geomagnetic  cutoffs,  at  r.uwn 
altitude,  as  reported  by  Shea  and  Smart  (1975). 
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LONGITUDE 

LATITUDE 

RIGIDITY 

LONGITUDE 

LATITUDE 

RIGIDITY 

80.00 

0.00 

0.02 

70.00 

0.00 

0.26 

80.00 

15.00 

0.04 

70.00 

15.00 

0.34 

80.00 

30.00 

0.06 

70.00 

30.00 

0.41 

80.00 

45.00 

0.09 

70.00 

45.00 

0.47 

80.00 

60.00 

0.09 

70.00 

60.00 

'  0.49 

80.00 

75.00 

0.10 

70.00 

75.00 

0.51 

80.00 

90.00 

0.10 

70.00 

90.00 

0.52 

80.00 

105.00 

0.11 

70.00 

105.00 

0.55 

80.00 

120.00 

0.11 

70.00 

120.00 

0.59 

80.00 

135.00 

0.11 

70.00 

135.00 

0.60 

80.00 

150.00 

0.08 

70.00 

150.00 

0.62 

80.00 

165.00 

0.06 

70.00 

165.00 

0.56 

80.00 

180.00 

0.04 

70.00 

180.00 

0.47 

80.00 

195.00 

0.00 

70.00 

195.00 

0.36 

80.00 

210.00 

0.00 

70.00 

210.00 

0.23 

80.00 

225.00 

0.00 

70.00 

225.00 

0.13 

80.00 

240.00 

0.00 

70.00 

240.00 

0.06 

80.00 

255.00 

0.00 

70.00 

255.00 

0.00 

80.00 

270.00 

0.00 

70.00 

270.00 

0.00 

80.00 

285.00 

0.00 

70.00 

285.00 

0.00 

80.00 

300.00 

0.00 

70.00 

300.00 

0.00 

80.00 

315.00 

0.00 

70.00 

315.00 

0.05 

80.00 

330.00 

0.00 

70.00 

330.00 

0.11 

80.00 

345.00 

0.00 

70.00 

345.00 

0.18 

75.00 

0.00 

0.10 

65.00 

0.00 

0.58 

75.00 

15.00 

0.14 

65.00 

15.00 

0.72 

75.00 

30.00 

0.18 

65.00 

30.00 

0.80 

75.00 

45.00 

0.20 

65.00 

45.00 

0.89 

75.00 

60.00 

0.23 

65.00 

60.00 

0.93 

75.00 

75.00 

0.25 

65.00 

75.00 

0.97 

75.00 

90.00 

0.25 

65.00 

90.00 

1.01 

75.00 

105.00 

0.26 

65.00 

105.00 

1  .03 

75.00 

120.00 

0.27 

65.00 

120.00 

1.12 

75.00 

135.00 

0.28 

65.00 

135.00 

1.19 

75.00 

150.00 

0.26 

65.00 

150.00 

1.20 

75.00 

165.00 

0.24 

65.00 

165.00 

1 .13 

75.00 

180.00 

0.20 

65.00 

180.00 

0.95 

75.00 

195.00 

0.14 

65.00 

195.00 

0.74 

75.00 

210.00 

0.09 

65.00 

210.00 

0.53 

75.00 

225.00 

0.03 

65.00 

225.00 

0.32 

75.00 

240.00 

0.00 

65.00 

240.00 

0.17 

75.00 

255.00 

0.00 

65.00 

255.00 

0.09 

75.00 

270.00 

0.00 

65.00 

270.00 

0.05 

75.00 

285.00 

0.00 

65.00 

285.00 

0.04 

75.00 

300.00 

0.00 

65.00 

300.00 

0.08 

75.00 

315.00 

0.00 

65.00 

315.00 

0.16 

75.00 

330.00 

0.02 

65.00 

330.00 

0.28 

75.00 

345.00 

0,07 

65.00 

345.00 

0.42 

LONGITUDE  LATITUDE  RIGIDITY  LONGITUDE  LATITUDE  RIGIDITY 


60.00  0.00  1.11* 

60.00  15.00  1.34 

60.00  30.00  1.46 

60.00  45.00  1.57 

60.00  60.00  1.61 

60.00  75.00  1.67 

60.00  90.00  1.73 

60.00  105.00  1.82 

60.00  120.00  1.95 

6C.00  135.00  2.05 

60.00  150.00  2.05 

60.00  165.00  1.99 

60.00  180.00  1.75 

60.00  195.00  1.40 

60.00  210.00  1.00 

60.00  225.00  0.65 

60.00  240.00  0.40 

60.00  255.00  0.22 

60.00  270.00  0.16 

60.00  285.00  0.14 

60.00  300.00  0.21 

60.00  315.00  0.38 


60.00 

330.00 

0.59 

60.00 

345.00 

C.86 

55.00 

0.00 

1.94 

55.00 

15.00 

2.28 

55.00 

30.00 

2.47 

55.00 

45.00 

2.61 

55.00 

60.00 

2.68 

55.00 

75.00 

2.78 

55.00 

90.00 

2.85 

55.00 

105.00 

2.92 

55.00 

120.00 

3.12 

55.00 

135.00 

3-31 

55.00 

150.00 

3.35 

55.00 

165.00 

3-15 

55.00 

180.00 

2.88 

55.00 

195.00 

2.22 

55.00 

210.00 

1.75 

55.00 

225 • 00 

1.23 

55.00 

240.00 

0.78 

55.00 

255.00 

0.50 

55.00 

270.00 

0.36 

55.00 

285.00 

0.36 

55.00 

300.00 

0.46 

55.00 

315.00 

0.75 

55.00 

330.00 

1.13 

55.00 

345.00 

1.59 

50.00 

0.00 

3-21 

50.00 

15.00 

3.54 

50.00 

30.00 

3.81 

50.00 

45.00 

3.97 

50.00 

60.00 

4.14 

50.00 

75.00 

4.27 

50.00 

90.00 

4.36 

50.00 

105.00 

4.37 

50.00 

120.00 

4.68 

50.00 

135.00 

4.93 

50.00 

150.00 

4.92 

50.00 

165.00 

4.67 

50.00 

180.00 

4.27 

50.00 

195.00 

3-38 

50.00 

210.00 

2.81 

50.00 

225.00 

2.03 

50.00 

240.00 

1.41 

50.00 

255.00 

0.95 

50.00 

270.00 

0.73 

50.00 

285.00 

0.69 

50.00 

300.00 

0.R9 

50.00 

315.00 

1.34 

50.00 

330-00 

1.98 

50.00 

345.00 

2.65 

li  r~t 

45.00  0.00  4.77 

45.00  15.00  5.12 

45.00  30.00  5.36 

45.00  45.00  5.51 

45.00  60.00  5.73 

45.00  75.00  5.90 

46.00  90.00  6.’1 

45,00  105.00  6.29 

45.00  120.00  6.j7 

45.00  135.00  6.86 

45.00  150.00  6.86 

45.00  165.00  6.33 

45.00  180.00  5.59 

45.00  195.00  4.85 

45.00  210.00  4.08 

45.00  225.00  3-16 

45.00  240.00  2.37 

45.00  255.00  1.74 

45  00  270.00  1.32 

45.00  285.00  1.22 

45.OO  300.00  1.49 

45.00  315.00  2.21 

45.00  330.00  3.16 

45.00  345.00  4.20 
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LONGITUDE 

LATITUDE 

RIGIDITY 

LONGITUDE 

LATITUDE 

RIGIDITY 

40.00 

0.00 

6.75 

30.00 

0,00 

1.30 

40.00 

15.00 

7.27 

30.00 

15.00 

1.71 

40.00 

30.00 

7.48 

30.00 

30.00 

2.  13 

40.00 

45.00 

7.70 

30.00 

45.00 

2.67 

40.00 

60.00 

8.19 

30.00 

60.00 

3.34 

40.00 

75.00 

8.73 

30.00 

75.00 

4.07 

40.00 

90.00 

9.14 

30.00 

90.00 

4.37 

40.00 

105.00 

9.29 

30.00 

105.00 

4.40 

40.00 

120.00 

9.49 

30.00 

120.00 

4.26 

40.00 

135.00 

9.89 

30.00 

135.00 

3.95 

40.00 

150.00 

9.74 

30.00 

150.00 

3.44 

40.00 

165.00 

8.95 

30.00 

165.00 

2.72 

40.00 

180.00 

7.86 

30.00 

180.00 

1.65 

40.00 

195.00 

6.46 

30.00 

195.00 

0.48 

40.00 

210.00 

5.41 

30.00 

210.00 

9.63 

40.00 

225.00 

4.55 

30.00 

225.00 

8.78 

40.00 

240.00 

3.61 

30.00 

240.00 

7.00 

40.00 

255.00 

2.76 

30.00 

255.00 

5.60 

40.00 

270.00 

2.07 

30.00 

270.00 

4.44 

40.00 

285.00 

1.93 

30.00 

285.00 

4.07 

40.00 

300.00 

2.42 

30.00 

300.00 

4.87 

40.00 

315.00 

3.41 

30.00 

315.00 

6.98 

40.00 

330.00 

4.82 

30.00 

330.00 

9.67 

40.00 

345.00 

5.92 

30.00 

345.00 

0.60 

35.00 

0.00 

9.54 

25.00 

0.00 

3. 10 

35.00 

15.00 

9.89 

25.00 

15.00 

3.64 

35.00 

30.00 

0.10 

25.00 

30.00 

4. 10 

35.00 

45.00 

0.53 

25.00 

45.00 

4.53 

35.00 

60.00 

1.15 

25.00 

60.00 

5.06 

35.00 

75.00 

1.44 

25.00 

75.00 

5.58 

35.00 

90.00 

1.52 

25.00 

90.00 

5.85 

35.00 

105.00 

1.71 

25.00 

105.00 

5.79 

35.00 

120.00 

1.93 

25.00 

120.00 

5.49 

35.00 

135.00 

2.04 

25.00 

135.00 

5.03 

35.00 

150.00 

1 .55 

25.00 

150.00 

4.44 

35.00 

165.00 

0.60 

25.00 

165.00 

3.76 

35.00 

180.00 

9.49 

25.00 

180.00 

3.07 

35.00 

195.00 

8.97 

25.00 

195.00 

2.43 

35.00 

210.00 

7.65 

25.00 

210.00 

1.81 

35.00 

225.00 

6.12 

25.00 

225.00 

0.98 

35.00 

240.00 

5.21 

25.00 

240.00 

9.74 

35.00 

255.00 

4.25 

25.00 

255.00 

7.89 

35.00 

270.00 

3.19 

25.00 

270.00 

6.08 

35.00 

285.00 

2.89 

25.00 

285.00 

5.44 

35.00 

300.00 

3.58 

25.00 

300.00 

6.56 

35.00 

315.00 

4.98 

25.00 

315.00 

9.65 

35.00 

330.00 

6.80 

25.00 

330.00 

1.47 

35.00 

345.00 

8.70 

25.00 

345.00 

2.46 
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LONGITUDE  LATITUDE  RIGIDITY 


LONGITUDE  LATITUDE  RIGIDITY 


20.00 

0.00 

4.11 

20.00 

15.00 

4.62 

20.00 

30.00 

5.09 

20.00 

45.00 

5.57 

20.00 

60.00 

6.12 

20.00 

75.00 

6.63 

20,00 

90.00 

6.87 

20.00 

105.00 

6.75 

20.00 

120.00 

6.37 

20.00 

135.00 

5.83 

20.00 

150.00 

5.20 

20.00 

165.00 

4.55 

20.00 

180.00 

3.93 

20.00 

195.00 

3.40 

20.00 

210.00 

2.91 

20.00 

225.00 

2.33 

20.00 

240.00 

1.47 

20.00 

255.00 

9.69 

20.00 

270.00 

7.82 

20.00 

285.00 

6.84 

20.00 

300.00 

8.65 

20.00 

315.00 

1 . 45 

20.00 

330.00 

2.76 

20.00 

345 .00 

3.54 

15.00 

0.00 

4.61 

15.00 

15.00 

5.14 

15.00 

30.00 

5.65 

15.00 

45.00 

6.17 

15.00 

60.00 

6.74 

15.00 

75.00 

7.25 

15.00 

90.00 

7.47 

15.00 

105.00 

7.34 

15.00 

120.00 

6.93 

15.00 

135.00 

6.37 

15.00 

150.00 

5.76 

15.00 

165.00 

5.17 

15.00 

180.00 

4.63 

15.00 

195.00 

4.16 

15.00 

210.00 

3.76 

15.00 

225.00 

3.30 

15.00 

240,00 

2.62 

15.00 

255.00 

1.27 

15.00 

270.00 

9.49 

15.00 

285.00 

8.46 

15.00 

300.00 

0.71 

15,00 

315.00 

2.52 

15.00 

330.00 

3.43 

15.00 

345.00 

4.07 

10.00 

0.00 

4.73 

10.00 

15.00 

5.26 

10.00 

30.00 

5.80 

10.00 

45 . 00 

6.34 

10.00 

bO .  00 

6.94 

10.00 

75.00 

7.44 

10.00 

90.00 

7.67 

10.00 

105.00 

7.56 

10.00 

120.00 

7.18 

10.00 

135.00 

6.65 

10.00 

150.00 

6.10 

10.00 

165.00 

5.61 

10.00 

180.00 

5.16 

10.00 

195.00 

4.75 

10.00 

210.00 

4.39 

10.00 

225.00 

4.00 

10.00 

240.00 

3.44 

10.00 

255.00 

2.55 

10.00 

270.00 

1.50 

10.00 

285.00 

1.07 

10.00 

300.00 

2.16 

10.00 

315.00 

3.09 

10.00 

330.00 

3.76 

10.00 

345.00 

4.25 

5.00 

0.00 

4.50 

5.00 

15.00 

4.99 

5.00 

30.00 

5.52 

5.00 

45.00 

6.10 

5.00 

60.00 

6.71 

5.00 

75.00 

7.22 

5.00 

90.00 

7.47 

5.00 

105.00 

7.42 

5.00 

120.00 

7.11 

5.00 

135.00 

6.66 

5.00 

150.00 

6.21 

5.00 

165.00 

5.83 

5.00 

180.00 

5.49 

5.00 

195.00 

5.15 

5.00 

210.00 

4.02 

5.00 

225.00 

4.47 

5.00 

240.00 

4.03 

5.00 

255.00 

3.41 

5.00 

270.00 

2.67 

5.00 

285.00 

2.40 

5.00 

300.00 

2.78 

5.00 

315.00 

3.36 

5.00 

330.00 

3.80 

5.00 

345.00 

4.11 
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LONGITUDE 

LATITUDE 

RIGIDITY 

LONGITUDE 

LATITUDE 

RIGIDITY 

0.00 

0.00 

3.94 

-10.00 

0.00 

2.11 

0.00 

15.00 

4.37 

-10.00 

15.00 

2.31 

0.00 

30.00 

4.87 

-10.00 

30.00 

2.71 

0,00 

45.00 

5.46 

-10.00 

45.00 

3.29 

0.00 

60.00 

6.10 

-10.00 

60.00 

3.89 

0.00 

75.00 

6.62 

-10.00 

75.00 

4.35 

0.00 

90.00 

6.90 

-10.00 

90.00 

4.68 

0.00 

105.00 

6.94 

-10.00 

105.00 

4.90 

0.00 

120.00 

6.73 

-10.00 

120.00 

4.91 

0.00 

135.00 

6.38 

-10.00 

135.00 

4.80 

0,00 

150.00 

6.05 

-10.00 

150.00 

4.74 

0.00 

165.00 

5.81 

-10.00 

165.00 

4.84 

0.00 

180.00 

5.59 

-10.00 

180.00 

4.94 

0.00 

195.00 

5.32 

-10.00 

195.00 

4.90 

0.00 

210.00 

5.03 

-10.00 

210,00 

4.76 

0.00 

225.00 

4.71 

-10.00 

225.00 

4.55 

0.00 

240.00 

4.33 

-10.00 

240.00 

4.29 

0.00 

255.00 

3.06 

-10.00 

255.00 

3.94 

0.00 

270.00 

3-32 

-10.00 

270.00 

3. 51 

0.00 

285.00 

2.95 

-10.00 

285.00 

3.10 

0.00 

300.00 

3.05 

-10.00 

300.00 

2.92 

0.00 

315.00 

3.38 

-10.00 

315.00 

2.82 

0.00 

330.00 

3.58 

-10.00 

330.00 

2.55 

0.00 

345.00 

3.69 

-10.00 

345.00 

2.20 

-5.00 

0.00 

3.13 

-15.00 

0.00 

0.75 

-5.00 

15.00 

3.45 

-15.00 

15.00 

0.91 

-5.00 

30.00 

3.91 

-15.00 

30.00 

1.23 

-5.00 

45.00 

4.50 

-15.00 

45.00 

1.69 

-5.00 

60.00 

5.14 

-15.00 

60.00 

2.18 

-5.00 

75.00 

5.65 

-15.00 

75.00 

2.70 

-5.00 

90.00 

5.97 

-15.00 

90.00 

3.00 

-5.00 

105.00 

6.10 

-15.00 

105.00 

3.25 

-5.00 

120.00 

6.00 

-15.00 

120.00 

3.25 

-5.00 

135.00 

5.77 

-15.00 

135.00 

3.32 

-5.00 

150.00 

5.58 

-15.00 

150.00 

3.48 

-5.00 

165.00 

5.50 

-15.00 

165.00 

3.78 

-5.00 

180.00 

5.42 

-15.00 

180.00 

4.09 

-5.00 

195.00 

5.25 

-15.00 

195.00 

4.24 

-5.00 

210.00 

5.02 

-15.00 

210.00 

4.24 

-5.00 

225.00 

4.74 

-15.00 

225.00 

4.15 

-5.00 

240.00 

4.41 

-15.00 

240.00 

3-97 

-5.00 

255.00 

4.01 

-15.00 

255.00 

3.69 

-5.00 

270.00 

3.53 

-15.00 

270.00 

3.30 

-5.00 

285.00 

3.14 

-15.00 

285.00 

2.89 

-5.00 

300.00 

3.08 

-15.00 

300.00 

2.59 

-5.00 

315.00 

3.19 

-15.00 

315.00 

2.31 

-5.00 

330.00 

3.15 

-15.00 

330.00 

1 .76 

-5.00 

345.00 

3.04 

-15.00 

345.00 

1.15 
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LONGITUDE  LATITUDE  RIGIDITY 


LONGITUDE  LATITUDE  RIGIDITY 


-20.00 

0.00 

9.21 

-20.00 

15.00 

9.06 

-20.00 

30.00 

9.29 

-20.00 

45.00 

9.70 

-20.00 

60.00 

0.21 

-20.00 

75.00 

0.25 

-20.00 

90.00 

0.41 

-20.00 

105.00 

0.65 

-20.00 

120.00 

0.84 

-20.00 

135.00 

0.62 

-20.00 

150.00 

0.74 

-20.00 

165.00 

1.93 

-20.00 

180.00 

2.78 

-20.00 

195.00 

3.19 

-20.00 

210.00 

3.44 

-20.00 

225.00 

3.51 

-20.00 

240.00 

3.46 

-20.00 

255.00 

3.28 

-20.00 

270.00 

2.94 

-20.00 

285.00 

2.53 

-20.00 

300.00 

2.14 

-20.00 

315.00 

1.65 

-20.00 

330.00 

0.73 

-20.00 

345.00 

9.78 

-25.00 

0.00 

7.55 

-25.00 

15.00 

7.50 

-25.00 

30.00 

7.72 

-25.00 

45.00 

7.88 

-25.00 

60.00 

7.86 

-25.00 

75.00 

7.37 

-25.00 

90.00 

7.08 

-25.00 

105.00 

7.24 

-25.00 

120.00 

7.43 

-25.00 

135.00 

7.71 

-25.00 

150.00 

8.47 

-25.00 

165.00 

9.55 

-25.00 

180.00 

0.22 

-25.00 

195.00 

1 .22 

-25.00 

210.00 

2.01 

-25.00 

225.00 

2.63 

-25.00 

240.00 

2.76 

-25.00 

255.00 

2.72 

-25.00 

270.00 

2.46 

-25.00 

285.00 

2.04 

300.00 

1.53 

315.00 

0.74 

330.00 

9.63 

S®5 

8.36 

-30.00 

0.00 

6.36 

-30.00 

15.00 

6.02 

-30.00 

30.00 

5.86 

-30.00 

45.00 

5.79 

-30.00 

60.00 

5.42 

-30.00 

75.00 

5.23 

-30.00 

90.00 

5.10 

-30.00 

105.00 

5.16 

-30.00 

120.00 

5.19 

-30.00 

135.00 

5.39 

-30.00 

150.00 

5.89 

-30.00 

165.00 

6.58 

-30.00 

180.00 

7.99 

-30.00 

195.00 

9.45 

-30.00 

210.00 

9.43 

-30.00 

225.00 

0.64 

-30.00 

240.00 

1.85 

-30.00 

255.00 

2.00 

-30.00 

270.00 

1.88 

-30.00 

285.00 

1.43 

-30.00 

300.00 
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